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Abstract 

Terrestrial planet atmospheres must be in long-term radiation bal- 
ance, with solar radiation absorbed matched by thermal radiation 
emitted. For hot moist atmospheres, however, there is an upper limit 
on the thermal emission which is decoupled from the surface temper- 
ature. If net absorbed solar radiation exceeds this limit the planet 
will heat uncontrollably, the so-called “runaway greenhouse”. Here 
we show that a runaway greenhouse induced steam atmosphere may 
be a stable state for a planet with the same amount of incident solar 
radiation as Earth has today, contrary to previous results. We have 
calculated the clear-sky radiation limits at line-by-line spectral resolu- 
tion for the first time. The thermal radiation limit is lower than previ- 
ously reported (282 W m~ 2 rather than 310 W m~ 2 ) and much more so- 
lar radiation would be absorbed (294 W m~ 2 rather than 222 W m -2 ). 
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Avoiding a runaway greenhouse under the present solar constant re- 
quires that the atmosphere is subsaturated with water, and that cloud 
albedo forcing exceeds cloud greenhouse forcing. Greenhouse warming 
could in theory trigger a runaway greenhouse but palaeoclimate com- 
parisons suggest that foreseeable increases in greenhouse gases will be 
insufficient to do this. 

Earth’s surface and troposphere are heated by the Sun. If more sunlight 
was absorbed, these would warm and, consequently, emit more thermal ra- 
diation to space. Whilst thermal radiation from the surface can escape the 
atmosphere directly, energy balance and temperate climate are maintained. 
Most such emission occurs through the “water vapour window”, the 8 to 
13 /jrn local minimum in water vapour absorption, coincident with the Wcin 
peak in surface thermal emission at 10 pm. 

Warming an Earth-like planet would make the atmosphere moist and 
hence optically thick across the thermal region, even in the water vapour 
window. Then, only the upper troposphere (not the surface or lower tro- 
posphere) would be able emit radiation directly to space. As sunlight pen- 
etrates the atmosphere better than thermal infrared radiation, the upper 
troposphere is heated from below and its thermal structure determined by 
moist convection. The water vapour mixing ratio increases with temperature, 
causing the moist adiabatic lapse rate to tend towards the saturation vapour 
pressure curve and the tropopause to acquire a fixed temperature-pressure 
structure. Thus the level from which effective thermal emission occurs tends 
to a fixed temperature and the flux a fixed value. This maximum in outgoing 
radiation — or “radiation limit” — means that surface warming no longer leads 
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to more thermal emission. If the net absorption of solar radiation exceeds 
this limit, then surface temperatures will increase in a “runaway greenhouse” , 
evaporating the entire ocean and sterilising the planet en route. 

This radiation limit was first found in a model by Simpson in 1927 
(ref. [1]), but seen as a paradox. The physics, as described above, was 
subsequently elucidated by Nakajima et al. [2], so we refer to the Simpson- 
Nakajima limit. In the intervening time, the notion of a runaway green- 
house was identified by reference to a separate stratospheric radiation limit, 
the Komabayashi-Ingersoll limit [3, 4] (which is unlikely to be reached in 
practice[2]). These various radiation limits have recently been reviewed[5]. 
Numerical calculations performed a few decades ago by several groups[6, 
7, 8, 9] converged on a consensus estimate of 310 W m -2 for the Simpson- 
Nakajima limit [9]. Subsequent work addressed some aspects of the problem 
(e.g. the hydrological cyclc[10], pedagogic treatment of thermal limits [11], 
early Mars[12]), but there has been no full revision of the classic calculations. 

Here, we present the most complete study of the runaway greenhouse for 
25 years, across the full range of temperatures, using modern input spectro- 
scopic data and a line-by-line treatment of the solar and thermal radiation 
(see the Methods section). We limit ourselves to clear-sky (cloud- free) cal- 
culations, which embody the first-order physics of the problem. Clouds both 
reflect solar radiation (making the runaway less likely) and enhance the green- 
house (making it more likely). Omitting them yields hard upper bounds on 
both solar absorption and thermal emission. This is a robust place to begin 
a re-evaluation of the problem, with thought experiments on how clouds will 
modify the results. 
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The runaway greenhouse has contemporary relevance. There has been 
high profile speculation that extreme anthropogenic global change could 
trigger it [13], but this is contrary to existing theory [2, 5] and numerical 
results[14]. It also sets the inner boundary of the circumstellar “habitable 
zone” [15, 16] in which the Kepler mission is presently identifying planetary 
candidates[17]. 


Pure water atmospheres 

We begin with the absolute end member case of a pure water atmosphere: 
neither any background gas nor any greenhouse gas other than water. With 
increasing surface temperature, evaporation from the ocean adds mass to 
the bottom of the atmosphere, so this is a good approximation for hot atmo- 
spheres in which the water vapour mixing ratio asymptotes to one. 

The optical depth of the atmosphere, r A (a function of wavelength), is 
measured downward from the top. Effective emission to space and atten- 
uation of sunlight occur where T\ ~ 1 (given Beer’s law; I\ = J 0 ex p -TA , 
where Jo and I\ are the incident and transmitted radiance), so plotting the 
altitude of T\ = 1 shows where emission to space and absorption of sunlight 
dominantly occur (Fig. 1) 

For thermal emission, T\ = 1 is either near the surface or not reached for 
low surface temperatures ( T s ), but “lifts off” towards high altitudes as the 
planet warms. Whilst T s < 1600 K, the temperature of T\ — 1 remains be- 
tween 250 and 300 K, so the top of atmosphere thermal spectrum is bounded 
by Planck functions for temperatures of 250 and 300 K, and is independent 
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of T s . Thus the Simpson-Nakajima radiation limit emerges as 282 Wm -2 
(Fig. 2), lower than previous estimates (e.g. 310 Wm -2 , ref [9]). When 
T s > 1600 K the upper atmosphere temperature gradient is sufficiently steep 
that the temperature reaches 400 K at T\ — 1 in the 4 /jrn water vapour 
window (Fig. lc). Hence a new peak in thermal radiation emerges, suffi- 
cient to permit a new stable climate with a steam atmosphere. No surface 
radiation escapes directly to space. Observation of this emission peak in a 
(exo)planetary atmosphere would indicate that the planet is in a runaway 
greenhouse state. 

Earth’s atmosphere is largely transparent to solar radiation. However, 
for water rich atmospheres, increasing temperatures are accompanied by in- 
creases in atmosphere pressure and water vapour absorption. The additional 
pressure increases the Rayleigh scattering optical depths at shorter wave- 
lengths, while near infrared water vapour vibration rotation bands increase 
the absorption optical depth at longer solar wavelengths. Both processes at- 
tenuate sunlight, very little of which reaches the surface, so surface albedo no 
longer affects the radiation budget. In the limiting case, for a pure water at- 
mosphere without clouds and the present solar flux, a maximum of 294 W m -2 
is absorbed, much higher than the previous estimate[9] of 222 W m -2 . 

Given a hot, moist and cloud free atmosphere, the net absorption of sun- 
light would slightly exceed the thermal radiation limit. This imples that 
a cloudless runaway greenhouse, steam atmosphere, would be stable under 
the present insolation. Earth today has a stable temperate climate (the 
requirements for which are discussed in the next section) implying a cli- 
mate bistability with respect to the runaway greenhouse (previously seen 
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in a grey atmosphere model [18]). Both the solar and thermal calculations 
represent upper bounds for a pure water atmosphere — clouds could reduce 
either, moving the bifurcation point. 

Previous work suggested that the thermal radiation limit does not depend 
on presence of non-condensible greenhouse gases[9, 2], This is not strictly 
correct. The radiation limit depends on the minimum absorption cross sec- 
tion in the 10 /mi water vapour window; any additional opacity here would 
raise the T\ = 1 surface to a higher altitude that radiates at a lower temper- 
ature. As purely theoretical tests, we set 1% each of our atmospheres to be 
carbon dioxide or ammonia, then 1% of both. The radiation limit decreased 
by 2 W m -2 for 1% C0 2 , 6 W m -2 for 1% NH 3 and 8 W m -2 for both (supple- 
mentary figure 7). Ammonia is one of the strongest absorbers around 10 /mi, 
so deeper reductions to the radiation limit seem unlikely. 

These results are sensitive to the absorption cross sections used. Firstly, 
using the most detailed spectral line-list for water (we use H1TEMP2010[19]) 
and correct Rayleigh scattering cross sections for water are essential. Using 
a less comprehensive line list (e.g. fflTRAN 2008), or Rayleigh scattering 
for air instead of water, gives erroneous results (supplementary figures 3 
6). Relative to previous results, our lower thermal emission and higher solar 
absorption are due to these absorption coefficient changes and increased spec- 
tral resolution. Secondly, the strength of the water vapour continuum (the 
smoothy varying absorption in the window regions) is very important. In 
the infrared, the continuum we use is weaker than indicated by the most 
recent data[20] (supplementary figure 9), so our estimates are conservative 
and the Simpson-Nakajima limit is likely somewhat lower than our estimate. 
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The uncertainty wedge associated with the continuum grows towards shorter 
wavelengths, and there are no measurements in the visible region. Our (or 
any other) solar calculations must be regarded as provisional until such mea- 
surements are made. 


Transition to a runaway greenhouse 

Given that a cloud-free steam atmosphere is a stable state at present so- 
lar constant, one should examine both how the stable temperate climate 
is maintained on Earth and the conditions which would lead to a runaway 
greenhouse. Hence we examined transitional atmospheres (up to 400 K) with 
the same mass of background gas as Earth and various greenhouse gas in- 
ventories (Table 1). 

At 280 K, the surface emits directly to space through the water vapour 
window (Figure 3). For surface temperatures above 310 K the temperature 
of the emitting level remains between 250 to 300 K, regardless of the surface 
temperature. If greenhouse gases other than water are more abundant, T\ = 1 
is higher in the absorption bands of these gases and less radiation is emitted 
overall. However, the relative magnitude of this effect decreases in hotter 
atmospheres with more water. In flux terms (Figure 4), for the end-member 
case of a saturated, cloud free atmosphere with contemporary surface albedo, 
the net absorbed solar radiation exceeds thermal emission in all scenarios 
except that with no greenhouse gases other than water, implying that a 
runaway greenhouse should occur. As this has manifestly not happened 
to Earth, we are led to the conclusion that a combination of atmospheric 
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subsaturation and an excess of cloud albedo forcing over cloud greenhouse 
forcing prevents a runaway greenhouse on Earth today. 

First, we relax the assumption of saturation. Our nominal relative hu- 
midity profile (see Methods) yields a stable climate only for the case of no 
additional greenhouse gases, although a marginally stable result was obtained 
for pre-industrial greenhouse gas concentrations. The assumed relative hu- 
midity profile is a source of uncertainty in 1-D models. With convection 
parameterisation in a 1-D model, Renno[10] found a non-linear transition 
from sub-saturation to saturation around 310 K, introducing an additional 
bistability in climate which is not found with climatological relative humidity 
profiles like ours[10]. In Earth’s tropics, columns of dry air prevent a local 
runaway greenhouse [21], Nonetheless, our results indicate that subsaturation 
alone is probably not sufficient to prevent a runaway greenhouse today (Fig. 

4). 

Thus, we turn our attention to clouds. Today, these give an albedo forc- 
ing of 50 Wm -2 and a greenhouse forcing of 26 Wm -2 (ref [22]). Climato- 
logical mean top of atmosphere fluxes are 239.4 W nr -2 net solar absorbed 
and 238.5 Wm~ 2 outgoing thermal[23] (the 0.9 Wm -2 discrepancy is caus- 
ing global warming). High, cold, clouds have a dominant greenhouse effect 
(though also reflect). Low cloud has a dominant albedo effect. At first 
approximation, we simulate this by increasing surface albedo; doubling our 
12% albedo approximately represents the low cloud forcing. Our subsat- 
urated preindustrial cases with present surface temperature and enhanced 
albedo (figure 4) have TOA solar and thermal fluxes of 263 and 266 W rri -2 , 
reasonably approximating climatology minus high cloud forcing. In all our 



scenarios (other than arbitrarily high C0 2 ), there is a stable climate under 
this assumption. Thus we can infer that the excess of cloud albedo over 
cloud greenhouse forcing, in combination with subsaturation, permits stable 
temperate climate on Earth. 

The clear sky fluxes are upper bounds: clouds could reduce either flux, 
making the runaway more or less likely. Previously [9] it was argued that 
cloud reflection would dominate over cloud greenhouse in an optically thick 
atmosphere. However, this misses a critical distinction based on the level 
of the clouds. The atmosphere is more transparent to sunlight than ther- 
mal radiation, so T\ (solar) = 1 is at lower altitude than T\ (thermal) = 1, so 
there are three categories of cloud effect based on these. Below T\ (solar) = 1, 
clouds will have negligible effect. Between t\ (solar) = 1 and T\ (thermal) = 1, 
albedo will dominate. Above T\ (thermal) = 1, cloud greenhouse will likely 
dominate (though reflection may dominate if the clouds are more than four 
times thicker than the present global mean[24]). The largest projected in- 
crease in the water vapour mixing ratio is in the upper atmosphere (supple- 
mentary figure 1), suggesting that high clouds would increase most (though 
this is speculative). For near future global warming, the current best esti- 
mates are for clouds to exert a positive forcing[25] (i.e. enhanced greenhouse 
dominating) . 

Steady state climates exist where the net flux (thermal minus solar) is 
zero; stable where the net flux increases with temperature and unstable where 
it decreases with temperature (Fig. 4). For small greenhouse gas inventories, 
the outgoing thermal flux overshoots the Simpson-Nakajima limit giving a 
“hump of stability: excess thermal emission will give a negative feedback, 
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restoring stable, temperature climate. Using the subsaturated, 25% albedo 
runs as a reference, the “hump of stability” is 24 Wm -2 for preindustrial, 
18 W m -2 for RCP 8.5 at 2100 and 8 W m -2 for extreme anthropogenic. With 
our arbitrarily high greenhouse gas scenario (30,000 ppmv C0 2 , this vanishes 
and there is no stable temperate climate. Greenhouse gases do not sim- 
ply warm the planet, but also lower or remove the energy barrier between 
temperate climate and a runaway greenhouse. 


Other times and other planets 

A runaway greenhouse has manifestly not occurred on post-ffadean Earth — 
it would have sterilised Earth (there is observer bias). Palaeoclimate gives 
us a sample of conditions where a runaway greenhouse did not occur, but 
cannot tell us the size of any safety margin. The “hothouse” climate of 
the Eocene is the most useful constraint for anthropogenic change. With 
solar constant 1% less than today and a few thousand ppmv C0 2 , the mean 
temperature was ~ 10 K warmer than today [26]. With C0 2 and temperature 
both higher then than we expect in the foreseeable future [27], this implies 
that an anthropogenic runaway greenhouse is unlikely. Deglaciaton from 
Neoproterozoic “snowball Earth” events likely required that ~10% of the 
atmosphere was carbon dioxide. The solar constant was 6% less than today, 
so net solar radiation absorbed would have been 12 W m -2 less and climate 
not yet bistable. By contrast, deglaciation from the next Snowball Earth 
might trigger a runaway. 

Venus probably experienced a runaway greenhouse in the past, evident 
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now in enrichment of D/H in its atmosphere [28, 7]. Previous work suggested 
that early Venus was close to the threshold for a runaway greenhouse[9] — our 
new lower radiation limit and enhanced solar absorption imply that, given 
the same amount of nitrogen in the atmosphere as Earth, Venus may not 
have had a habitable period. However, if early Venus had at least as much 
nitrogen in its atmosphere as it does now this would have had a protective 
effect (Figure 5). More nitrogen gives more Rayleigh scattering, decreasing 
absorbed solar radiation in a transitional atmosphere. 

As the solar constant increases with time, Earth’s future is analogous to 
Venus’s past. We expect a runaway greenhouse on Earth 1.5 billion years 
hence if water is the only greenhouse gas, or sooner if there are others. Earth’s 
atmospheric nitrogen inventory has likely changed with time [29]. Any future 
decrease would lessen the protective effect of Rayleigh scattering and hasten 
a runaway greenhouse (somewhat by contrast to previous arguments [30]). 

Our pure water calculations were aimed at hot atmospheres, but should 
also apply to “water-worlds”, analogous to a warm version of Jupiter’s moon 
Europa. In the absence of Rayleigh scattering from background gas, plane- 
tary albedo would be lower than Earth. Under Earth’s insolation, without 
clouds and with 12% surface albedo (the average for Earth, about twice that 
of seawater), there is no stable temperate climate and runaway greenhouse 
would always ensue. Arbitrarily increasing surface albedo to 25% (a proxy 
for low cloud reflection) gives a marginally stable state at 275 K. With a 
mean surface temperature this low, ice albedo feedback would likely lead to 
low latitude glaciation (Snowball Earth-like). A transient warming sufficient 
to melt the ice would likely cause a transition directly to a runaway green- 
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house. Unless mediated by other atmospheric constituents or clouds, there 
would be no stable temperate climate state. 

The runaway greenhouse sets a hard limit for the inner edge of the cir- 
cumstellar “habitable zone” [15, 16]. This classic definition neglects multi- 
ple climate equilibria so may be misleading. Three major stable climate 
states exist at the same solar constant: Snowball Earth (at least transiently 
reduced habitability), temperate (habitable) or runaway greenhouse steam 
atmosphere (sterilizing) [18]. Changes to clouds today may be sufficient to 
transition from temperate to either other state. Determining surface tem- 
perature requires knowledge of the atmospheric state and history: it is not 
possible to determine a habitability a priori from incident stellar radiation. 

The thermal radiation limit depends weakly on the mass of the planet 
(Supplementary figure 8). Everything else being equal, a Mars size planet 
would be more susceptible to a runaway greenhouse and a so-called “super- 
Earth” less [11]. 

Revisiting the classic planetary sciences problem of the runaway green- 
house with modern modelling tools, we have shown that the thermal radiation 
limit is lower and that more solar radiation is absorbed. The runaway green- 
house may be much easier to initiate than previously thought. A renewed 
modelling effort is needed, addressing both Earth and planetary science appli- 
cations. We have begun this process with a single column, clear-sky model, 
which has allowed us to advance the core radiative transfer aspect of the 
problem. Reference calculations are available as Supplementary Information 
to permit future model testing. Our work should be followed with cloudy 
column models then, ultimately, general circulation models (to address the 


12 



cloud and relative humidity distributions). The latter represents a “grand 
challenge” in climate modelling, for which present-generation models may be 
insufficient: there are difficulties associated with radiative transfer, clouds 
and dynamics (with a major component being condensible), and no empiri- 
cal comparison cases. 
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Methods 


Atmospheric structure 

We prescribe the surface temperature and atmospheric structure (moist adi- 
abatic, with bottom of atmosphere temperature equal to surface tempera- 
ture) in a single global mean column [9]. We separately consider pure water 
atmospheres (which is the limit for extensive evaporation of the ocean and 
allows us to easily include the non-ideal behaviour of steam) and transitional 
atmospheres (up to a surface temperature of 400 K where the ideal gas law 
remains a good assumption). More information is given in the supplementary 
information, along with profiles (supplementary figure 1). 

Relative humidity 

The standard assumption for one- dimensional models is a linear decrease in 
relative humidity with pressure, from the surface to the tropopause[31, 32], 
Mean subsaturation on Earth today is an area-weighted construction of areas 
of saturated and unsaturated air. Once water vapour is a major atmospheric 
constituent (a few percent or more), regions of sub saturation would induce 
huge pressure gradients, leading to rapid mixing from saturated regions and 
permitting further evaporation. Hence we assume that air below a level with 
some threshold water vapour mixing ratio (we use 5%) is saturated, with a 
linear decrease in relative humidity with pressure above. 
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Spectral data 

We use line data for water from the HITEMP2010 spectral database[19], 
which provides line data, including weak lines, in the range 0-30,000 cm -1 
(>0.33 /irn; Supplementary figure 3) and absorption cross sections for the far 
UV (> 50, 000 cm -1 , < 0.2 /irn). The line data is converted to cross sec- 
tions separately for each atmospheric profile examined using the LBLABC 
program, written by David Crisp. There is no spectral data available in 
the mid UV (30000-50000 cm -1 , 0.2-0.33 /rm). This corresponds to 6% of 
incoming solar radiation, so we will underestimate net absorption. It is of 
utmost importance to use the most detailed line database available for water 
(supplementary figures 3-4). For thermal emission, HITRAN2004 and later 
perform similarly to HITEMP2010 up to 625 K whereas HITRAN2k allows 
more outgoing radiation through the 4 and 10 /irn windows. For solar radia- 
tion, all HITRAN datasets give less absorption than HITEMP2010. At high 
temperatures (thick atmospheres) the largest differences are seen at shorter 
wavelengths as net absorption depends on the ratio of absorption to Rayleigh 
scattering cross sections (the latter depending on the reciprocal of the fourth 
power of wavelength). For greenhouse gases other than water, we use line 
data from HITRAN2008[33]. 

Continuum absorption 

The water vapour continuum is smoothy varying absorption which cannot be 
explained by near-centre contributions of known spectral lines. It is primar- 
ily responsible for setting the absorption cross section in the water vapour 
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window, which in turn largely determines the radiation limit. There has 
been interesting theoretical work on the source of this recently [34], but it 
will be some time until this provides a mature parameterisation at the range 
of temperature, pressure and mixing ratio conditions that we encounter. Our 
approach is to use a more standard implementation, which is implemented 
with y-factors (ref. [35] and Clough, personal communication, 2000) and per- 
forms similarly in tests to the empirical MT-CKD 2.4 continuum. We then 
test our computed absorption cross sections against the best available data in 
the thermal infrared [20] (supplementary figure 9). Our implementation some- 
what underestimates the strength of the continuum, so the radiation limit 
may be lower than our estimate and thermal exit from the runaway would 
require a warmer temperature. Whilst the temperature-pressure range of 
the data is limited, we take comfort in that the thermal radiation limit is 
determined in the cool upper atmosphere, which is within observable ranges. 

Rayleigh scattering 

Previous work[9] used Rayleigh scattering cross sections for air throughout, 
as data for water was not available (J. Kasting, personal communication). 
Detailed refractive indices for water are now available [36] , so we derive new 
constants for Rayleigh scattering (supplementary table 1). Water is a weaker 
Rayleigh scatterer than air, so using Rayleigh scattering cross sections for 
air in place of water leads to erroneously low solar absorption, as too much 
sunlight is scattered (supplementary figure 6). 
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Radiative transfer calculation 


We use the SMART code, written by David Crisp[37], for our radiative trans- 
fer calculations. This code works at line-by-line resolution, but uses a spectral 
mapping algorithm to treat different wavenumber regions with similar opti- 
cal properties together, giving significant savings in computational cost. We 
evaluate the radiative transfer in the range 50-100,000 cm^ 1 (0.1-200 /an) as 
a combined solar and thermal calculation. Our solar source is spectrally re- 
solved and we average the flux from zenith angles of 7.5°, 22.5°, 37.5°, 52.5°, 
67.5°and 82.5°. In a few cases, models for one of these zenith angles did not 
run successfully due to singularities in the matrix inversion in DISORT [38], 
which SMART uses for multiple scattering. In these cases, an appropriate 
average was made of the remaining five and the effect on the final results 
is trivial. Spectra displayed are reduced via a simulated slit function. Our 
reference-standard radiative transfer model output is available as supplemen- 
tary information to facilitate future model testing. 
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Scenario 

C0 2 (ppmv) 

CH 4 (ppmv) 

N 2 0 (ppbv) 

Baseline 

0 

0 

0 

Preindustrial 

287 

0.806 

275 

RCP 8.5, year 2100 

936 

3.75 

435 

Extreme anthropogenic 

3000 

10 

500 

Arbitrarily high 

30,000 

10 

500 


Table 1: Greenhouse gas inventory scenarios used in transitional at- 
mospheres. Preindustrial [39] and RCP (representative concentration path- 
way) 8.5 for year 2100 (ref [40]) are standard. Extreme anthropogenic as- 
sumes that the large fossil fuel reservoirs (coal and some methane clathrate) 
are burned rapidly — i.e. it is extreme, but certainly possible. Note that con- 
centrations are relative to a standard, 10 5 Pa atmosphere and the mass of 
each greenhouse gas is fixed as temperature and surface pressure increase, so 
actual mixing ratios decrease as temperature increases. 


25 



Surface solar radiance down (W m * |im' 


a 





).2 1 2 
Wavelength Jim) 





T = 400K 


300 

200 

100 

0 

400 

300 

200 

100 




T = 300K 


500 1000 1500 

Atmospheric temperautre (K) 


,.i.i A A At.. 

.2 1 10 30 

Wavelength Jim) 



2 10 30 

Wavelength (gm) 


Figure 1: Spectra and effective absorption/emission levels in a pure 
water atmosphere T s is surface temperature, a. Spectrum of downward 
solar flux at the surface. No line data is available for the grey shaded area. b. 
Altitude at which optical depth is unity. Solid line is absorption optical depth 
and dashed line is Rayleigh scattering optical depth. Background shading 
is atmospheric temperature, c. Spectrum of outgoing thermal flux at the 
top of the atmosphere. Solid lines are the black body flux at the surface 
temperature. Dotted lines are the black body flux for 280 K and 400 K for 
comparison. See also supplementary figure 2 
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Figure 2: Top of atmosphere fluxes from a pure water atmosphere. 

Colours identify different surface albedos: green for 12% and purple for 25%. 
Steady state climates are found when the net outgoing flux is zero; stable 
steady states are where the flux is increasing with increasing temperature as 
it passes through zero and unstable steady states are where it decreases with 
increasing temperature as it passes through zero. 
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Figure 3: Spectra of thermal emission level and outgoing thermal 
radiation for transitional atmospheres Black is baseline case, red is with 
5000 ppmv C0 2 . a. Altitude where optical depth is unity. Most outgoing 
thermal radiation is emitted from this level. Background colour is atmo- 
spheric temperature, b. Top of atmosphere emission spectra. Solid grey 
lines are the surface emission for each case. Dotted grey lines are for Planck 
functions of 220 K, 250 K, 280 K, 310 K, 340 K and 370 K up to the surface 
temperature. The broad 5 to 8 /irn absorption feature is from water and the 
smaller absorption features at 4.2 /irn and 15/rnn are from carbon dioxide. 
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Figure 4: Top of atmosphere fluxes from transitional atmospheres. 

a-c are with a surface albedo of 12%. d-f are with a surface albedo of 25%. 
Colours identify the greenhouse gas inventories from Table 1 (black is base- 
line, bine is preindustrial, green is RCP 8.5 at 2100, red is extreme anthro- 
pogenic and purple is arbitrarily high). Solid lines are atmospheres saturated 
with water vapour throughout whereas dashed lines are subsaturated where 
the saturation mixing ratio is less than 5%. Steady state climates exist where 
the net outgoing flux is zero, stable where flux increases with temperature, 
unstable where flux decreases with temperature. 
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Figure 5: Top of atmosphere fluxes from an ideal gas atmosphere 
with a varying amount of background gas (nitrogen). All runs for 
baseline (no non-condensible greenhouse gas) with surface albedo of 12%. 
Colours are: blue for 0.1 bar, turquoise for 0.33 bar, black for lbar, orange 
for 3 bar and red for lObar. 
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